NACA TN 3018 


AUTHOR'S PERSONAL COPY f 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 3018 


|A THEORETICAL STUDY OF THE EFFECT OF FORWARD SPEED ON THE 


ту 


">p 
S 


FREE -SPACE SOUND-PRESSURE FIELD AROUND PROPELLERS 
By I, E, Garrick and C. E. Watkins 


Langley Aeronautical Laboratory 
Langley Field, Và, 


Š і тор, REFERENCE 
` ` BRARY COPY ун» 
JM i 986 yor то BETES 
з ЫН 
Washington 


October 1953 


1F 


LE ШШШ 


3 1176 01433 4305 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3018 


A THEORETICAL STUDY OF THE EFFECT OF FORWARD SPEED ON THE 
FREE-SPACE SOUND-PRESSURE FIELD AROUND PROPELLERS 


By I. E. Garrick and C. E. Watkins 
SUMMARY 


The.sound-pressure field of & rotating propeller in forward flight 
in free space is snalyzed by replacing the normal-pressure distribution 
over the propeller &ssociated with thrust and torque by a distribution 
of acoustic pressure doublets acting at the propeller disk and subject 
to uniform rectilinear motion. The basic element used to synthesize the 
field is the pressure field of a concentrated force moving uniformly at 
subsonic speeds, for which an expression generalizing one of Horace Lamb's 
for the fixed concentrated force is given. This result is presented both 
for the moving and for the fixed observer. The strength of the doublet 
distribution is related to the thrust and torque distribution in a con- 
venient but approximate way. The sound field is expressed by integration 
over the propeller disk, and also by integration over an effective ring, 
and is given both for the near pressure field and, in a simpler form, 
for the far field. Known results for the zero-forward-speed case present 
themselves in the special case of Mach number M= 0. Some illustrative 
examples are calculated and discussed. 


INTRODUCTION 


The rotating propeller is the source of an intense sound=pressure 
field which can be associated with the periodic reactions on the medium 


arising from the distribution of pressure rotating along with the blades. 


This pressure distribution consists in part of & distribution due to 
thickness of the blades, whose resultant force in subsonic potential flow 
is zero, and in part a distribution due to angle of attack and camber 

of the blades whose integrated effect includes the induced drag and corre- 
sponds almost wholly to the thrust and torque distribution over the blade. 
Another source of propeller noise may be associated with flow separation 
and with friction or shear due to the boundary layer; both effects lead 
to vorticity shed into the wake and hence the designation vortex noise. 
The vortex noise and the noise due to thickness (where wave drag is not 

& large factor) are, however, for actual propellers normally of а consid- 
erably smaller magnitude than the rotational sound due to torque and 
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thrust; hence only the latter effect will be considered in the present 
work. | 


А large number of investigators have studied various phases of the 
determination of the sound or noise field of rotating propellers. In 
addition to the references cited, a bibliography is included of repre- 
sentative work on this subject. A simplification that has frequently 
been made in propeller-noise investigations is to limit the consider- 
ations to the static or standing propeller. The work of Gutin in 1936 
(ref. 1) represents а development-of this type that makes possible a 
satisfactory prediction of the amplitude of sound pressure due to thrust 
and torque of a propeller rotating on a stand in still air. Although 
Gutin's theory is applicable to the near oscillating-pressure field of 
the propeller, his results in reference l are limited to the determi- 
nation of the fundamental and the first few harmonics at a distance far 
from the propeller, that is, several diameters away. The determination 
of the near pressure field, however, has been of concern both from struc- 
tural and physiological considerations. Hubbard and Regier (ref. 2) 
extended the application of Gutin's work to describe the oscillating 
pressure field and to determine the amplitude of noise at points near 
the propeller, in some cases within a blade chord length from the tip. 
They investigated analytically the effect of several of the parameters 
that enter in the theory and also gave comparisons with experiment which 
were quite satisfactory. 


The existing theoretical work has found useful application for 
static conditions and for conditions of low forward speed, for example, 
near take-off. For conditions of high forward speeds, however, many 
pertinent questions have arisen as to the possible effects of the forward 
speed on the oscillating pressure field of the propeller. А few investi- 
gators have examined phases of this problem; Bryan, Hart, Shirokov, 
Blokhintzev, and particularly Kuessner and Billings (see bibliography ) 
may be mentioned, but most of this work seems incomplete or difficult to 
apply. It appeared desirable, especially for applications, to reexamine 
the theoretical problem Гог the propeller at forward speed во as to have 
it arise as а straightforward generalization of existing work for the 
zero-forward-speed case. The purpose of the present paper is therefore 
to extend the theory of Gutin (ref. 1) and the work of Hubbard and 
Regier (ref. 2) so as to include the effect of subsonic forward speed 
of the propeller on the near and the far oscillating sound-pressure 
field caused by torque and thrust of the propeller. 


This paper includes the following material: (а) The sound-pressure 
field associated with a uniformly moving concentrated force is given with 
details of the explicit development contained in the appendix. This 
result serves as а basic element for synthesizing the total pressure 
field. It is expressed both for the case of the observer or field point 
considered to be moving along rectilinearly and uniformly with the рго- 
peller in free space and for the case of the observer considered fixed 
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and the propeller in uniform flight. (b) The disturbance forces associ- 
&ted with the pressure distribution &cting on the medium in the plane of 
the propeller disk are presented in the manner similar to that of Gutin. 
(c) The sound-pressure field resulting from the combined use of (а) and 
(b) is shown in the form requiring integrations over the propeller disk 
and in the simpler form making use of an effective propeller radius. 
Approximations valid for the far field, which reduce to Gutin's results 
for the case of zero forward speed, are also given. (d) Some numerical 
examples are calculated and discussed. (e) Remarks on the usefulness 
and limitations of the analysis. 


SYMBOLS 
A(r) chordwise distribution of thrust acting on a radial element 
of a propeller blade 
B number of propeller blades 
b(r) width of propeller blades 
c velocity of sound 
Cp power coefficient, Р/ра?ро 
Ст thrust coefficient, Т/оп2р" 
D propeller diameter 
F(r) chordwise distribution of forces perpendicular to the thrust 


of a propeller blade and giving rise to the torque 


Fy, F, components of force vector F 


х? су 
к= = 
c 
тв Bessel function of first kind with index mB 
M Mach number, У/с 
m order of harmonic 


їі propeller rotational speed, rps 


4 NACA TN 3018 


Р power 
p pressure 
pressure magnitude 
Po pressure due to torque 
Рт pressure due to thrust 
Pus root-mean-square pressure 
torque E 
R length of propeller blades 
Re effective length of propeller blades 
r radius to a blade element 
r, § polar coordinates in yz-plane 


S = (x - ху)? + B? (y ~ yi) + (z - 2d 


бо = x2 + e82 000 | й u 


8 = (х = xi) + (y у) + (® Е 21)? 


во = хе + y? 


Sos 9 polar coordinates in xy-plane 
T thrust 

t time 

У forward velocity 


x, y, 5 Cartesian coordinates 
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BR blade angle 
T propeller efficiency 
p fluid density 

M(x - ху) + 8 
g = —————ə— ' h 

82 

р 
т = — 

го 
Ta period, 2x/BQ 
Q angular velocity 
O) frequency of mth harmonic, mBQ 
© fundamental frequency, ВО 


ANALYSIS 
Extension of a Formula in Lamb's "Hydrodynamics" for the 
Sound-Pressure Field of a Fixed Concentrated Force 
to That of a Moving Force 
On the basis of acoustic considerations of the classical hydrodynamic 


equations, Lamb (ref. 3) gives the pressure at any field point x, y, z 


associated with an external periodic force Ae lot acting in the x-direc- 


tion and concentrated at the location хі, уу, 21 ав 


- 8. 
PIE 8 


where 


в = (x - m+ (у - у) + (z = зі)" (2) 
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The concentrated force may be regarded as stemming from the limit of a 
distribution of an increasing pressure difference over a decreasing 
area, whose product in the limit of zero area is equal to the force. 


For an arbitrary concentrated force of components F, Fy» F, the result 
generalizes as Lamb indicates to 
8 5 8 
зайво; 18) a (ts), ә F(t - с) (3) 
ha | Ox 8 ду B oz S 


This formul& shows tbat the pressure field obtained from the concentrated 


force has the character of a doublet or a dipolel; in its periodic form 
it has played a central role in the development of the noise and pressure 
fields of propellers due to the propeller torque and thrust. However, 
the formula refers to a concentrated force and to coordinates fixed in 
space. To treat the case involving effects of forward speed, it is expe- 
dient to obtain the required extension of the preceding formula. 


The extension of equation (3) required for the case of a uniformly 
moving concentrated force is given in equation (4), which follows. 
Details for its derivation are supplied in the appendix. Let the con- 
centrated force have components F,, Fy» Р.) let it be moving uniformly 


with velocity V іп the direction of the positive x-axis; &nd let the 
coordinate system also move uniformly with the same velocity. Then, 

with x, у, ® now denoting the field point referred to this coordinate 
system and with ху, yy, 21 denoting the coordinates of the concentrated 


moving force я a 


G С С 
je aja Ке а Ве - 2) miss 
Ця | Ox 5 ду 5 Oz S 
where К. ds - T" 
S = V(x - хіў + Sa (y - у)? + (2 - "92 
"m - xy) + 8 (5) 
Be 
в = үз - м 
lsince the acoustical pressure is p= -o & the dipole is an accel- 


ot’ 
eration doublet or the time rate of change of а flow doublet. 
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Figure 1 illustrates the geometric meaming ої the quantities в, 5, 
and о. Тһе force located at  0(х),у1,21) and the field point at 


Q(x,y,z) are both in uniform motion with velocity У in the x-direc- 
tion. The distance OQ = s. The influence at Q at time t stems 

from the action of the force when it was at position P, where P is 
obtained from the relation oP/0Q = M, the forwerd-speed Mach number. 

The distance PQ gives then directly the phase radius c. Те perpen- 
dicular dropped from О onto PQ determines QR, which is equal to the 
amplitude radius S. For Ме O both S and с reduce to the ordinary 
radius s. Comparison of equations (5) and (4) shows that the effect of 
the forward speed leads to replacement of s by 5 in the amplitude and 
by of in the phase. 


The force variation of main concern herein is that of a harmonically 
periodic force F having components Ру, Fy, F, varying in time as eiot, 


Equation (4) becomes 


à à à JT 


> (дн ФЕ = (б) 


hat х У ду 2 Oz 


where К = а/с. Equation (6) is the basic relation to be made use of іп 
the subsequent analysis. 


Disturbance Forces in the Propeller Plane 


Consider as does Gutin (ref. 1) а propeller disk oriented so that 
the axis passing through the center of the disk coincides with the x-axis 
and let the propeller be assumed to rotate, as in figure 2, in the 
yz-plane (x = O) with positive values of x corresponding to points 
ahead of (and negative values of x corresponding to points behind) the 
propeller disk. The propeller is considered to move uniformly with 
velocity V in the positive x-direction. Points in the propeller disk 
are designated by (0534524); Or, in polar coordinates (as shown in 


fig. 2), by 


Уу = г cos 8 


(T) 


21 = r sim 8 


For definiteness, let the propeller be rotating counterclockwise as seen 
ру an observer looking into the propeller toward the slipstream. 
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Each element of the propeller is acted on by the surface pressure 
distribution and this distribution may be resolved into a thrust force 
in the direction of the axis of rotation (the x-axis) and into а force 
associated with the torque which acts about the axis of rotation opposing 
the rotation. Equal and opposite reaction forces to these are exerted on 
the medium. The points of application of these forces are imagined to 
act in a single plane designated as the plane of rotation or as the pro- 
peller disk. | | Т | 


The reaction of the surface pressure distribution of the rotating 
propeller on the medium at any instant is to be replaced by fixed 
periodic forces acting at the propeller disk. The propeller disk itself 
will be considered covered with the necessary singularities in the pres- 
sure or acoustic radiators of proper strength and harmonic content to 
correspond to this normal-pressure distribution. These singularities 
will be seen to be acceleration sources (acoustic radiator of zero order) 
for those forces which act symmetrically on both sides of the blade, that 
is, whose net force over each element of the blade is zero, and to be 
acceleration doublets (acoustic radiator of order one) for those forces 
acting antisymmetrically, that ls, whose net force over each element of 
the blade corresponds to the difference in pressure over both sides of 
the blade. The doublet distribution is that which is needed to repre- 
sent the thrust and torque distribution, in particular that part of the 
thrust and torque distribution associated with pressures acting normal 
to the blade surface. This part is practically all of it, except that 
arising from skin friction. (The effect of blade thickness may be taken 
into account by introducing either flow sources or acceleration sources: 
the flow sources would lead to the sound-pressure field due to thickness 
by the classical "piston" effect of the moving blade on the flow field; 
the acceleration sources would deal directly with the assumed or known 
contribution to the pressure distribution over the blade due to thickness.) 


Consider an element of the propeller at distance r from the axis; 
let dr be its radial length and b its width measured in the projec- 
tion onto the plane of rotation. Let the forces acting on the propeller 
element on each blade be A(r)dr in the axial flight direction and 
F(r)dr in the direction opposed to the direction of rotation. Equal 
reaction forces acting opposite to the flight direction and in the direc- 
tion of rotation, respectively, are exerted on the medium. These quanti- 
ties are related to thrust T and torque 0 by the relations 


= = BA(r) 


(8a) 
= = BrF(r) 


2F 
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or 


A(r)àr 


а 

il 
E 
E пу 


(8b) 


rF(r)dr 


D 
T 
tJ 
o 


where B ів the number of blades and R the length of each blade. 


The periodic impulse or reaction experienced at any element of the 
disk may be expanded in a Fourier series. For simplicity let the element 
considered first be located on the radial line Ө = О (the y-axis) and 
afterwards be located arbitrarily. To start with, assume that the forces 
are uniformly distributed over the projected width b of the blade 
element, that is, that the distribution of pressure difference over the 
blade chord is rectangular. Then in the area element r dr dO of the 


plane of rotation, the forces A(r)dr = 22 and Р(г)аг = 32 


medium during the time interval іп which this element is eclipsed by the 
projection of the propeller element. If the overlapping of the element 


starts at t= 0, it will end at + = т = 2 and the overlapping of the 


act on the 


element will start again by the next blade at + = To = = where Q is 


the angular velocity of the propeller. The rectangular-type forces 
experienced at the element of the disk located at 0 = O by its periodic 
eclipse by the blade may be developed in a Fourier series: 


г dé 


A(r) dr for Oc t<T 
FJ ($) = b 
O for T< t< To 
x со 
F. (t) = Ао + > An cos(mBQOt - Em) (9) 
m=1. 
and similarly 
x F(r) 2-89 ar for O<t<T 
Fo (t) = b 
O for T< t< Tg 
со 
Fa (t) = Bo + > Bp cos (mBQt - Ти? (10) 
m=1 
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where 
2 г T 
= — A(r) — sin mx — dr dé 
Am mx ( ) b To 
2 г T 
‚ = — Fir) — sin mt — dr ds 11 
Ва (г) 2 = (12) 
_ T _ шВОт 


The constants As and Во which correspond to the instantaneous average 


thrust and torque over the blade element and to the associated momentum 
shed into the slipstream do not give rise to sound and need not be 
expressed. Те phase angles є, and Ту; which are small for the lower 


harmonics, are needed to preserve generality of discussion since A, cos €m 
corresponds to the coefficients of pure cosine terms and А, sin €, to 


coefficients of pure sine terms in the Fourier series. For the assumed 
rectangular distribution, є and Ta may be observed to be equal to 
zero if the origin of time is chosen at the overlapping of the center 
line of the blade by the y-axis since then only a pure cosine series 
suffices to express the distribution. 


Ihe general formulation of interest, where the disk element r dr дө 
is located at an angle 8, may be expressed directly from these results, 
for, on a second area element г аг ад shifted with respect to the first 
by the angle @ in the rotational direction, there act periodic forces of 
the same magnitude, but retarded by the time в/о. The corresponding 
Fourier developments are | 


F(t) = Fy gi - 2) = Ag + > Am соз(шВоё - шВӨ - em) 
(12) 
Fo(t) = то" Í - =) = Bo к. Вл cos(mBQt - шВе - n4) 
The quantity 
элт — mBOT — mBb 
m a2 on (13) 
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is small Гог the lower harmonics, especially Гог blade elements not near 
the hub. Blade elements near the hub, where r is small, are eliminated 
from consideration since they contribute very little to the air forces. 


One may then replace in equation (11) sin 1 by = so that, noting 
о о 


the relations in equations (8a), 


Ag = A(r)àr ao = = ar ав 


(14 ) 
B. 140 


This approximation, which weights the harmonic content of all the har- 
monics equally, may be observed to correspond to the case of the thrust 
and torque distributed over a zero blade width (that is, to the mathe- 
matical pulse sometimes termed the Dirac delta function). This approxi- 


mation becomes relatively less valid when T exceeds about і 
О 
(sin = O. TOT; E = 0.785) or when hm — exceeds about unity; that is, 
O 


when the order of the harmonic exceeds = or, roughly, 1/solidity. 


Moreover, the assumption of a uniform rectangular distribution of 
the forces across the blade has been made for convenience, and other 
distributions may be treated if desired. Some of the possible errors in 
the assumption leading to equation (14) are discussed briefly by Gutin 
in reference 1 and are shown to be generally negligible for the lower 
harmonics. Regier and Hubbard (ref. 4) also discuss this assumption in 
an illuminating manner. Figure 3, which is taken from reference h, 
shows the relative harmonic content of different assumed distributions 
of the same total load for: (a) the mathematical sharp pulse, (b) a 


triangular “hat" pulse ae 0.03, (c) a rectangular pulse 8 0.03, 
2nr 2xr 


and, (d) a rectangular pulse z^ 0.06. Blade widths of actual pro- 


pellers tend to be between cases (c) and (d) for the most effective parts 
of the propeller disk. The sharp-pulse assumption generally tends to 
overestimate the magnitude of the higher harmonics. When the proper 
distribution is known, appropriate correction factors for the required 
harmonic may, of course, be applied to the magnitude of the results given 
by the pulse solution. 
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The formulation in equation (12) need not be limited to the one in 
which the propeller force distribution is uniform throughout the pro- 
peller cycle. If interest should be attached. to & nonuniform distribu- 
tion, &s occurs for the propeller yawed or pitched with respect to the 
flight path or for one experiencing interference effects, it is readily 
possible to &llow for these effects by permitting the amplitude of the 
distribution to become а function of 9. 


The Sound-Pressure Field 


The general expressions for the resolution of the forces associated 
with the thrust and torque оп a radial blade element (for example, 
eqs. (12)) may be put in the usual convenient complex form (whose real 
part may correspond to the formulation of interest): 


EN NES Es a o (et=mB8-cm) 
(15) 
Fo = Bo + > В 
m=1 


where w= mB and where for the special case of the rectangular pulse 
distribution the coefficients are given as in equation (11). 


With the use of the approximations for A, and B, given in equa- 


tion (14+), the periodic forces acting on the medium at the element г dr 48 
of the disk can be expressed for any given harmonic in the forms 


_ Lal i(wt-mBe) 
_ 2149 1 (wt-mBe ) 

Fy m = = аз singe dr dé (16) 
. 1149 1 (аб -пів8 ) 


In the derivation ої equation (16), the phase angles єп and 1, which 


can generally be neglected for the lower harmonics, have been put equal 
to zero, a value which corresponds in the case considered (as has been 
mentioned) to the choice of Ө = О to correspond to the overlapping of 
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the center line of one of the blades at the time ^t = О. The index m 
indicates the order of the harmonic considered, Ав discussed in the 
preceding section, appropriate factors depending on the harmonic num- 

ber m, functions of the radius, and, for unsymmetrical loading, func- 
tions also of Ө may (if known) be applied to the terms in equation (16). 


Pressure relations involving integrations over the propeller disk.- 
The pressure at any point in the free-space field produced by these com- 


ponents of the uniformly moving periodic forces F,, Fy, F, is given in 


equation (6). With the use of the components of force given in equa- 
tion (16) and by integration over the propeller disk, the total oscil- 
lating pressure p for any given harmonic m (the index m will be 
hereinafter dropped) is obtained from equation (6) as a sum of the pres- 
sure đue to thrust рт and that due to torque Ро 88 follows: 


р = Ppt Po (17) 
where 
R nex -iko 
Pp = -— L QT S 2 e ас ae (18) 
ha^ Јо Jo & 8 
R 2x -iko 
Pg = - 1 | 1 49 ,i(wt-mBe) 1 д e —Ó (19) 
un-Yo Уо ra гд 8 
and where, in the expression for Po» the operation == has been used 


` for convenience to replace the equivalent operation sin ° = - cos ө È. 


The indicated differentiations with respect to x and Ө сап be carried 
out by use of the expressions for c and S given in equations (5) 


ОЕ КО. er Ao Die ж 
Әх Sg S 


(20) 


-ikc -iko 2 
Loe eg eB SD ig NN. 
roe s S 5 g2 
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Hence, 


iot pR penx -іка | 
„= dT Q-imBO 9 "7 ey Eu am 5" do (21) 


and 


iot pR рол -iko 2 
р. = 8 ЕІ = 99 .-imbe e "(ik |, P. (y sin 0 - z сов 0)àr do 
Q ще ovo * б NS ge M" т | 


(22) 


Another somewhat simpler expression Гог Po can be obtained from equae- 
tion (19) by integration by parts with respect to 6 


imB dr 48 (23) 


eint ГА 2л 1 aq -impo е-іко 
r 


Pressure relations involving effective ring approximation. - Appreci- 
able simplification may be achieved for calculation purposes by making 


use of the approximations inherent in the assumption of an effective pro- 
peller radius БК. so that the integration with respect to r is avoided. 


Equations (21) and (23) then reduce to 


iwt pen -i (mB9+ko ) 
рр = e J (ER. ikx | x e де (ац) 
hx? Jo 82 Bes g2 S 


and 


lot per -i (mB6+ko 
= / З imB атікшвонко) do (25) 
O 


Ree 8 


where, in both S and с, the points J4 and 21 һауе the values 


Уу = Re соз Ө, 2) = Ва sin Ө. In effect, equations (24) and (25) imply 


that the propeller disk has been replaced by an annuler ring in which 
the entire thrust and torque are concentrated. The effective radius of 
this ring varies somewhat, as calculation shows, with the load distri- 
bution and with the order of the harmonic. Deming (ref. 5) has shown by 
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calculation in special cases for the static propeller that the ring 
approximation is a reasonably good one. An effective radius of the order 
of O.8R is considered reasonable for adoption in initial calculations. 


The magnitude of the root-mean-square pressure Diana is of interest, 


Since most sound-recording instruments are calibrated in terms of it. The 
contribution of the harmonic of order m to the root-mean-square pres- 
sure is 


г 


2 
pum = |Рт + Pol 


or 
Pong = 125002 + вам? (26) 
8x 
where 
ex | ту k x mB ав 
к= j == cos(mBO + ko) + т (M+ 2)- @ 28. [sinus + xo) бо 
o 182 p? S Re? S 


2л | 
B= |. -IX sin(mBO + ко) + т S e =) - q P. |cos(mBO + kc) 22 


When rotational symmetry exists as for the condition of symmetrical 
loading, it is convenient in making the numerical calculations to let 
the field point be in the xy-plane, an arrangement which can always be 
attained by suitable choice of the line Ө = O. 


Pressure relations. for the far field.- A further simplification in 


the results can be reached if the distance s from the propeller disk 
(0,1, 21) to the field point (x,y,O) is large. Then from 


e 


S 2х2 (у зу) + 24° 


there is obtained 


в = sç - — (27) 
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where So = хе + y2 (the distance from the field point to the center 
of the propeller disk). It also follows from 


62 = x? +—B2(y - уз)? + 82272 
that 


(28) 


where 
Зо = Vx? + Bey? 


If terms of order 1/s? are neglected in comparison with terms of 
order 1/8, the sound pressure due to thrust is obtained from equa- 
tion (24) by use of equations (5) and (28): 


-İ— (So +Mx 
рр = Т ° Р вон) EU + a аа Lands ae (29) 
Т hx? So ве 50/У о | 
From the known integral relation 
est : 
O 
it follows that 
-1— [S +Мх 
ете YB. 2 ome) зм + e (31) 
Т 2x8, ge So So 
Similarly, from equation (25), 
г E eit mB+l -1— (SotMx ) (Se) (32) 
Q RS? 2x99 Sg 
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Observe that the argument of the Bessel function may be replaced by 


ы. = ВМ. оь — (55) 
O O 


QR 
where Meot = =, the Mach number corresponding to the rotational speed 


at the effective radius. (The far-field approximation and the effective- 
radius approximation made use of in arriving at eqs. (51) and (32) need 
not, of course, be made simultaneously, since the required integration 
of the Bessel functions with respect to the radius may be carried out, 

at least numerically, without difficulty. ) 


Introduce the fundamental frequency denoted by ау = BQ, so that 


The pressure magnitude for any harmonic m is then given by 


Pl = [Pp + ва| 
moy XX Вс kyRe 
= ———— — — - J. 
PEA "(м + р ен "UE: М 5 (5%) 


This result шау be compared with that for the case of zero forward speed 
(M = O) given in reference 1: 


ШОУ; 
2nCSo 


[P| = aae) (35) 


It can be noted that the forward-speed Mach number affects the 


torque term containing 9 only in the replacing of в, = ух? + y2 (the 
distance from the field point to the propeller hub) by the smaller dis- 


tance Зо = үх2 + p2y2. This substitution occurs both in the argument 
of the Bessel function and in the outside factor 1/S,. The thrust term 


containing T is more strongly affected, since, in addition to this 
change, the field effect ahead of and behind the propeller disk is influ- 


enced by the term M+ = 


and also by an increase because of the factor 1/82. 


corresponding to а backward shift by = = -М 
о 
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APPLICATION ТО А SPECIFIC PROPELLER 


In order to. give some indication of the effect of forward-speed 
Mach number on the sound pressure of a propeller, calculations based on 
equation (26) for the near field and on equation (34) for the far field 
have been made. For this purpose, a two-blade research propeller having 
а lO-foot diameter and operating under the various conditions summarized 
in table I was chosen, As may be noted from this table, the propeller 
is assumed to operate at constant power; that is, the-power coeffi- 
cient Ср and the torque Q are held constant. As the forward-speed — 
Mach number is changed, the blade angle Вв and consequently the | 


thrust T are changed according to the propeller charts of reference 6 
so as to make the assumed conditions consistent with actual test oper- 
ating conditions. Only the fundamental, the first harmonic m = 1, is 
considered in these examples and the value chosen for the effective 
radius Ве is О.8К or 4-feet. It should be pointed out that the sound 


pressure computed from the data of table I would be obtained in pounds 
per square foot, since these data are given in English units. These 
values have been converted to dynes per square centimeter in the results 
given in the figures by multiplication by the approximate conversion 
factor 480. 


Calculations for the near field.- Calculations of the root-mean- 


square pressure (Ps = ХА мірі) based. on equation (26) are made for 
various values of x in the range from -0.5D to O. 5D along the line 

y = 0.6D (6 feet), that is, along a line 1/5 of the propeller radius from 
the tip and extending a distance of 1 radius behind to 1 radius ahead of 
the plane of rotation. These results are shown plotted as a function of 
x/D for several Mach numbers in figure 4. It can be noted in this fig- 
ure that, for each Mach number, two peaks generally appear, one ahead of 
and one behind the plane of rotation. For the вер ої conditions under 
consideration (see table I), the highest peak amplitudes of pressure 
occur just behind the propeller plane at values of x/D in the range 
from -0.15 to -0.075. Аз the Mach number is increased from О to 0.4, 

the peak amplitudes decrease in magnitude, but as the Mach number is 
increased from 0.4 to 0.9 this trend is reversed and the peak amplitudes 
increase in magnitude, the peak for М = 0.9 being about 1.4 times that 
for Ме 0. Also, ав M increases from O to 0.9, the point at which 
the highest peak pressure occurs moves somewhat nearer the propeller 
plane. Thus, the generally severe sound-pressure conditions at take-off 
(M = O) tend to be alleviated in flight at the lower Mach numbers but may 
be reached again and even exceeded at the higher flight Mach numbers. 

The calculations for low Mach numbers presented in figure 4, it should 
be noted, are in substantial agreement with both calculated and measured 
results of reference 2. The high peak pressures obtained for the highest 
Mach numbers indicate that, with propeller-driven airplanes, the sound 
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pressures generated near the tips of the propellers are of significance 
and of possible concern with regard to both structural considerations 
and passenger comfort. 


Calculations for the far field.- Calculations of the root-mean- 
square pressures based on equation (34) for various values of x along 
the line y - 2D (20 feet) are shown plotted in figure 5. The trend with 
regard to Mach number in this figure is about the same as was noted in 
figure 4 for y = 0.6D, but the relative effect of the forward Mach num- 
ber appears greater at a distance. However, because of the greater dis- 
tences in cds 2 than in figure 4, the peak Бирез Luss are considerably 
less. 


Directional characteristics of the sound field.- For some purposes, 
especially with regard to calculations for the far field, it is desirable 
to consider the character of the sound pressure in terms of polar coordi- 
nates. For this purpose, substitutions can be made in equation (3L) 

(and eq. (26)) as follows: 


о cos Ф у = so sin Ф So = Ух? + y2 (36) 


The quantities x/S, and y/S, in equation (34) (and eq. (26)) 
may then be written as 


x _ сов Ф 
So V1 - M?sin^o 
(5T) 
D e sin 6 
8 


Vi - МевіпеФ 


With these substitutions, equation (34) becomes 


cos Ф Be KRe sin Ф 
-======='\- Q ——|x | — 
Vi - = этне) алВе (1 - M2sin2ç 

(38) 


Calculations based on equation (58) for the Mach numbers of O and 
0.8 for a constant value of в, = 2D are shown plotted as dashed curves 
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in figure 6. The solid curves in this figure represent results of calcu- 
lations along the line у = 2D obtained by replacing sg by = 

sin : 
with y = 2D in equation (50) for the various Mach numbers considered 
in figures 4 and 5. А comparison of results of calculations along the 
line у = 2D with those along the circle зо = 2D for M= 0 and 


М = 0.8 indicates that; for high Mach numbers (of the order of M = 0.8), 
the peak pressures calculated along a line y = Constant- are about the 
same as those calculated along the circle with radius equal to the con- 
stant value of у. Observe the second pressure peak which has developed 
in the forward location at M = 0.9. 


Separate components due to torque and to thrust of the sound field.- 
In order to give some indication as to the nature and proportion of 


sound pressure associated with each of the quantities, thrust T and 
torque Q, the root-mesn-square pressures associated with each of these 
quantities for М = O and М = 0.8 are plotted for у = 0.60 (6 feet) in 
figure T. These plots show, as do equations (26) and (34), that the root- 
mean-square pressures associated with the propeller torque are symmetri- 
cally distributed with respect to the plane of the propeller for all Mach 
numbers, whereas the root-mean-square pressures associated with the 

thrust are symmetrical with respect to this plane only Гог M= O. For 
the particular propeller and operating conditions under consideration, 

the amplitudes of pressure associated with thrust for low Mach numbers 

are higher than those associated with torque, риб for high Mach numbers, A 
the opposite is true. Тп the interpretation of figure T, it should be 

recalled that the results depend on the assumed operating conditions and 

that the torque Q is the same (2,680 ft-1b) in both parts of the figure; 
whereas the thrust is 1,850 pounds at М = O and 310 pounds at М = 0.8. 

Plots of this type can be used to obtain the sound pressure for various 

thrust and torque coefficients for а given propeller, since these coef- 

ficients appear as factors in equations (26) and (34) and hence can be 
normalized. 


It is of particular interest to note in the plot Гог М = 0.8 that 
the peak pressures associated with both thrust and torque are-consider- 
ably greater than the peak pressure associated with the sum of these 
quantities. This result indicates that the phase relationship between 
the two components is, in this case, such that each has a canceling 
effect on the other. Also in this case, the pressure associated with 
thrust may have its greatest value ahead of the plane of the propeller; 
whereas the total pressure has its greatest value behind the plane of 
the propeller. | 


The discussion of. the figures based on the specific examples illu- 
strates the fact that the results do depend markedly on the assumed 
operating conditions. Moreover, only the fundamental harmonic (m = 1) : 
for the specific two-blade. propeller has been illustrated. Trend studies 
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on effects of higher harmonics, number of blades, amd different орег- 
ating conditions would be of considerable interest. Some preliminary 
calculations on the sound pressures associated with the fundamental of 
а six-blade propeller have shown а greater relative effect of the 
forward-speed Mach number. 


CONCLUDING REMARKS 


Expressions have been given for the sound-pressure. field due to the 
distribution of thrust and torque for any given harmonic of a rotating 
propeller in uniform subsonic flight. Тһе general expressions (eqs. (17) 
to (25)) involve integrations over the propeller disk; approximate 
expressions for the near field (eqs. (2h) to (26)) and for the far field 
(eqs. (31), (32), and (34)) involve integrations over a ring with an 
effective radius. The numerical examples have illustrated some free- 
space sound-pressure results for the fundamental of a specific two-blade 
propeller under various operating conditions at various forward speeds. 


It is pertinent to remark again on some of the limitations of the 
analysis. The analysis presented utilizes the torque and the thrust 
distributions, which may be given empirically or theoretically, in such 
& way as to require that they &rise purely from pressures acting normal 
to the surface of the blades. Empirical values of the thrust and torque 
include & contribution, generally small, due to skin friction and to 
separated flow. Some caution is then needed in the use of the results. 
For example, the empirical torque term is somewhat larger than the torque 
due purely to normal pressures and, hence, the associated sound result 
due to torque may be overestimated slightly; similarly, the empirical 
thrust term шву slightly underestimate the sound due to thrust. At high 
tip speeds there may be significant contributions to the noise due to 
wave drag associated with the thickness. These contributions are mainly 
taken into account by the effects of the wave drag on the torque. (Other 
sound effects of the thickness should be separately calculated and 
included but these effects have not been explicitly presented herein. ) 


Another way of looking at the approximation, and perhaps a generally 
more convenient one for the study of trends, is to consider that the 
assumed torque and thrust distributions are actual theoretical ones, 
obtained if necessary by adjustment of the blades of the propeller іп а 
potential flow; and hence the sound-pressure field is that corresponding 
to the chosen thrust and torque, but in potential flow. Although the 
sound-pressure field and the aerodynamic velocity field have been con- 
sidered as separate, it is of interest and of significance that the same 
concepts leading to the calculation of the sound field (made use of in 
the form of the acceleration potential) can lead to the linearized aero- 
dynamics of the propeller in compressible flow, including the represen- 
tation of the vorticity left behind in the wake as a result of spanwise. 
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variation of loading. The theoretical induced drag and the theoretical 
wave drag (which is a form of acoustic loss) are inherently included in ` 
the representation employed. 


It may be worthy of repetition that the pressure formulas as given 
tend to overestimate the contribution of the higher harmonics and that 
appropriate factors based on chordwise loading can be devised and may 
be required. Moreover, because of the theoretical (as well as the 
empirical) change in aerodynamic loading along the blade radially with 
increasing tip Mach number and forward-speed Mach number, the appropriate 
factors will change. In addition, the effective radius will be altered 
with the loading. Thus, although 0.8 blade radius may be suitable as an 
effective radius for normal loading, & smaller value for the effective 
radius may be more suitable for conditions where unloading of the tip 
occurs. = 


A calculation study of trends under different assumed conditions, 
the effect of unsymmetrical loading, or of dual propellers, effects of 
obstacles or boundaries on the free-space results, and experimental 
confirmations for the in-flight propeller are interesting matters for 
further investigation. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., August 31, 1953. 
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APPENDIX 


THE SOUND-PRESSURE FIELD ASSOCIATED WITH 
A UNIFORMLY MOVING FORCE 
Derivation of Equations (4) and (6) 
Lamb (ref. 3, p. 502) derives from hydrodynamical principles the 
differential equation satisfied by the pressure condensation s* in the 


acoustic field associated with arbitrarily varying fixed forces acting 
on the medium 


2 2 2 
с126% = > ° ELE 
32 дуг За 2 52 
_ 1/9х, дү, Z 
x. nm z) (A1) 


where Lamb's X, Y, Z are actually the extraneous forces per unit volume 
divided by the density p. (In effect, Iamb is dealing here with the 


&cceleration potential.) In terms of the perturbation pressure р = pc^g* 
the equation is 


сор = div F 
oF, oF, oF, 
s d NE ae VN 


where F is the arbitrary force per unit volume having components 
Fy, Fy, F,. Lamb shows that if the periodic force Ае10 is imagined 
concentrated on an infinitely small space at (хо зУ1»21,) and to be in the 


direction of x, the pressure at the field point (x,y,z) (the distance s 
from the location of the force) is given by 


where 
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so that the concentrated force is equivalent to a double source or 
acceleration doublet whose axis is in the direction of the force. For 


a general harmonically periodic force F having components Fx, Fy, Fz 


varying in time as et” 
„= LF. ç) Š 
peg (куа 
1 à Š 9 \ етіке 


where k= >. 


For an arbitrary time-dependent concentrated force F(t) located 
at (6,1,6), the pressure at field point (x,y,z) сап be expressed as 


p S 
doe EE 
ila F(t -2) à (6-3) ә (Е-е) 
то” "uae sea "mS иш лыш И. 


where the differentiations affect the result only through the variable в. 


The extension to a moving concentrated force may, within the frame- 
work of small-perturbation theory, be made in several ways - for example, 
as mentioned in reference 7. One formal procedure utilized by Kussner, 
for example, makes use of the invariance properties of the wave equation 
and utilizes combined Galilean and Lorentz transformations. The following 
direct procedure given briefly in reference 7 is believed to be of 
intrinsic interest; as it represents a simplification in Prandtl's pro- 
cedure (ref. 8) for the case of moving constant source distributions, a 
procedure which consists of scheduling & suecession of fixed sources in 
a, path to act consecutively one after the other во as to represent in 
effect the desired source moving along its path. Let the arbitrary con- 
centrated force act only as an impulse during an infinitesimal interval 
at time t= T. The impulse may be written as 


F(t)8(t - T) (A6) 


tF 
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where the impulse function б(т) = О for т # О and is characterized 
for т = О as having unit area with respect to T. The useful property 
of the impulse function of "sifting or selecting" a value of a function 
is exhibited by the following relation (see, for example, ref. 9, p. 61): 


Í f(r)5(c - т)дт 


- СО 


F f(c - т)% (т)ат 


== 


f(c) (АТ) 


Let а succession of such impulses act, one following the other, in а path, 
points along which are given in space-fixed coordinates by 


Е = Е(т) п = n(T) б = t(T) 


The effect at time t of all such impulses which act before the time t 
is then given from Lamb's result (eq. (A5)) as 


p= baw [BO (вот sje: (A8) 


where, as has been defined, there is & nonzero contribution to the 
integral only for values of T defined by the characteristic relation 


t-r = Sadie - кт] + [y - alr]? [s - £I? (A9) 


Which expresses the distance between the source point and the field point 
in terms of the time of tr&vel of the outgoing waves. 


The integral in equation (A8) corresponds to a summation of temporary 
fixed sources. To represent the case of uniform rectilinear motion with 
velocity V in the positive x-direction, let the sources be located on 
the €-axis and flow consecutively one after the other at the positions 


Ë = Vr T = 0 C= 0 (A10) 


so that for T = -æ the source is located at Ë = -о and for т = 0 
the source is located at the origin. It follows that the distance 
between source point and field point is 


8 = Vx = yr)? + у2 + z2 (А11) 
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Replace in equation (АЗ) the variable т by 9 where 
8 
t-T-<=°e (A12) 


By the sifting property of the 8-function (eq. (AT)) the integral will 
have a value only for 0 = O and equation (A8) becomes 


TUE 
rm atv] FC Ja (A15) 


where the quantity within the bracket is now to be evaluated. Equa- 
tions (All) amd (A12) give a quadratic equation for T which, on choice 
of the solution that leads to T< t and with 9 = 0, gives 


r= 3(¢- 4-8) TUS 
where 


8 = V(x - vt)? + 2(y2 + 22) 
В = Vl - м 


id 
c 


Equations (A11) and (A12) also give 


1 дт с 
S 48 es - V(x - Vr) 


where replacing s by its characteristic value c(t - т) from equa- 
tion (A9) and replacing т by its value given in equation (Al}) yields, 
for Ө = O, I 


E ucl 
do 8 (A15) 


where S is defined in equation (All), Equation (А15) becomes 


(A16) 
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This result is referred to space-fixed coordinates; it is finally 
desired also to express this result in terms of a field point (хо›Уо› 20) 


of а coordinate system moving uniformly along with the source located at 
ху, Уу; 21 so that one may put 


x - Vt = хо - Хі y = yo - Yı 2 = Zo - 2) 


After this substitution is made, the zero subscripts may be dropped to 
yield as the end result corresponding to equation (A5) for the source 
and the coordinate system in uniform rectilinear motion: 


11 
= i _\__°/ 
a= div 7 (A17) 
where 
S = (x - хт)? + S2 (y - y1)* + (z 2 2)? 
M(x - X1) + 8 
Gg = —————— 
g2 
в «1 - м 


Equation (AlT) corresponds to the result given in equation (4) of 


the analysis. For the periodic harmonic force F varying as elot, ії 
may be expressed іп the form 
= -ikc 
p = = (F. ç) Š (A18) 
list 


which is equivalent to equation (6) of the analysis. Equations (A17) 
and (А18) are the sought-for generalizations of equations (A5) and (A) 
of the appendix, for the case where the disturbance and the field point 
are in uniform rectilinear motion. Comparison of these equations shows 
at once that s is replaced by S in the amplitude and by o in the 
phase. A geometric interpretation of these quantities is shown in fig- 
ure 1 and discussed in the analysis following equation (5). 
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Remarks On the Case of the Moving Disturbance 
and the Fixed Observer 


A few remarks are in order on the significance of equation (A16). 
In this equation, the field point РИ ів given іп space-fixed coordi- 
nates while the disturbance force (or propeller) is in the in-flight 
condition moving with velocity V in the positive x-direction. Its 
location is given by ху = Vt, yy = О, 2) = O so that at t-0 it is 


located at the origin. It may be observed that the distance between 
disturbance and field point S is numerically the same whether given 
by the relation in equation (A16) or (А17). Hence, the pressure magni- 
tude at any observer location for а sound-radiating element of the 
in-flight propeller is the same as that for the observer moving along 
with the propeller, provided the proper instantaneous distance between 
observer and propeller is used. However, there will be a difference in 
the frequency perceived by the observer. This frequency will be that 
for the case of the uniformly moving observer modified by the Doppler 
effect; thus, the frequency об each harmonic is modified by с (е + Vr); 


where У, is the component of the propeller.forward speed (with proper 


sign: minus, for &pproaching, pius for receding from, the observer ) 

in the acoustical direction from the observer to, say, for the far 
approximation, the hub of the propeller. The acoustical direction is 
not quite that from the observer to the propeller location; it actually 
points from the observer to the location of the propeller when the sound 
which reaches the observer was emitted (direction ОР rather than 90 
in fig. 1). The Doppler frequenegy factor e /( с + Vr) is given geomet- 


rically in figure 1 by the ratio of 5 = =. 
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TABLE І 
SEA-LEVEL OPERATING CONDITIONS FOR A 10-FOOT-DIAMETER 


TWO-BLADE PROPELLER (DATA FROM REFERENCE 6) 


E = 815 hp; Cp = 0.10; Q = 2,680 lb-ft; k = 0.29686 | 
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Figure 1.- Geometric representation of в, 5, and с. 
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Figure 2.- Propeller disk and coordinate system. 
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Figure 5.- Effect of impulse shape on the relative amplitudes 


of the harmonics. 
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Figure 4.- Distribution of the root-mean-square pressures for the 
fundamental (m = 1) of а two-blade lO-foot-diameter propeller at 
several forward-speed Mach numbers. Cp = 0.10; k = 0.29686; 


y = 6 feet. (Operating conditions of propeller given in table I.) 
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Figure 5.- Distribution ої the root-mean-square pressures for the 
fundamental (m = 1) of a two-blade 10-foot-diameter propeller at 
several forward-speed Mach numbers. Cp = 0.10; k = 0. 29686; 


y = 20 feet. (Operating conditions of propeller given in table I. ) 


38 NACA TN 3018 


100 дума/ст (а) Ма O; T=I85O b. 100 дюн/ст? (b) M = O.2, T=IGOO b. 
90° 
l2 
te) o | 
I8O O° І8 o? 
00 desem? (с) M = O.4, T=875 lb. ЮО dynes/em® (4) M = O.6; T=580O b. 


Flight direction ——Ə 
90° 


конуса? (f) M = OS; T=2O5 b. 


— — Line y*20 feet | м 


---- Circle s,=2O feet 


Figure 6.- Polar diagrams of the root-mean-square pressures calculated 
along a line у = 20 feet. and along а circle s, = 20 feet fora 


two-blade, 10-foot-diameter propeller at several forward-speed Mach ` 
numbers. Cp = 0.10; m = 1; k = 0.29686. (Operating conditions of 


propeller given in table I.) 
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T = 510 11. 


(b) M = 0.8; 


М = 0; T = 1850 lb. 
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Figure 7.- Distribution of the root-mean-square total pressures and the root-mean-square pres- 


(Operating con- 


› 10-foot-diameter propeller at forward- 


Cp = 0.10; m= l; k = 0.29686; y = 6 feet. 


ditions of propeller given in table I.) 


gures associated with thrust and torque for а two-blade 


speed Mach numbers of 0 and 0.8. 
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